The spectral form factor near the Ehrenfest-time 
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We calculate the Ehrenfest-time dependence of the leading quantum correction to the spectral form 
factor of a ballistic chaotic cavity using periodic orbit theory. For the case of broken time-reversal 
symmetry, our result differs from that previously obtained using field-theoretic methods [Tian and 
Larkin, Phys. Rev. B 70, 035305 (2004)]. The discrepancy shows that short-time regularization 
procedures dramatically affect physics near the Ehrenfest-time. 
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From a statistical point of view, the spectra of a quan- 
tum particle confined to a cavity with disorder or con- 
fined to a cavity with ballistic and chaotic classical dy- 
namics are remarkably similar 0. In both cases, the 
spectral statistics are independent of system specifics and 
follow the predictions of random matrix theory (RMT) 
on energy scales comparable to the mean spacing be- 
tween energy levels A or, equivalently, on time scales 
comparable to the Heisenberg time tn = 2nh/A 0, E|, 
provided that the latter is much larger than the time 
T crg needed for ergodic exploration of the phase space. 
Calculations of the spectral statistics have been the 
theoretical vehicle through which the profound corre- 
spondence between RMT, impurity diagrammatic per- 
turbation theory, periodic-orbit theory, and the field- 
theoretic zero-dimensional a model has been demon- 
strated HSHH0]. 

In the last decade, it has been understood that chaotic 
quantum systems are characterized by a time scale in- 
termediate between r org and th that does not have its 
counterpart in diffusive systems. This time scale is the 
'Ehrenfest time' r E |1 II El El The Ehrenfest time 
is the time required for two classical trajectories initially 
a quantum distance (wavelength) apart to diverge and 
reach a classical separation (system size) . It is expressed 
in terms of the Lyapunov exponent A of the correspond- 
ing classical dynamics as te = A -1 \n(c 2 /fi), where c 2 is 
a classical (action) scale. The existence of the Ehrenfest 
time does not affect the universality of the spectral statis- 
tics. However, it is responsible for differences between 
otherwise universal properties of chaotic and disordered 
quantum systems for energies ~ 7t/rR , which is well in- 
side the universal range if h/c 2 -> [H El dill El Cl- 
in this article we address the spectral statistics at this 
energy scale. 

Most of the literature on spectral statistics consid- 
ers the spectral form factor K(t), which is the Fourier- 
transform of the two-point correlation function of the 
level density p(e), 
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Here the brackets (. . .) c denote the connected average 
obtained by varying the center energy e and/or other pa- 
rameters in the system. For T crg <C t < in, the RMT 
predicts that K(f) is dominated by a perturbative ex- 
pansion in t/tn 
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where f3 = 1 (2) in the presence (absence) of time- reversal 
symmetry and 



SK 2 (t) 



0. 
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For times t > t-a the perturbative expansion in t/tn 
breaks down, and K(t) is governed by non-perturbative 
contributions [l6|. 

The presence of the Ehrenfest time does not affect 
the leading contribution to K(t), but it does impact 
SK . The leading te dependence of 5K% in the pertur- 
bative regime, which already occurs to order (t/tn) 2 , 
was first considered by Aleiner and Larkin [T^ . using 
a field-theoretic approach. Recently, Tian and Larkin 
used the field-theoretic approach to calculate the leading 
TE-dependence of SK 2 0, which appears to order (t/tu) 3 
and causes the perturbative quantum correction 5K% no 
longer to be strictly zero. Their result is 
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5K 2 {t) = - 



2irht 2 i 



9(t - 2r E ) + . . . , (4) 



[6(f - 3r E ) - 0(t - 4r E )] + . . . , (5) 



where 9(x) = 1 if x > and 9{x) = otherwise, 
0(x) = Jq dx'6{x') = x6(x), and the dots ... refer to 
terms of higher order in t/tn that were not considered 
in the calculation. (Tian and Larkin also considered the 
TE-dependence of non-perturbative contributions to the 
spectral form factor, but these will not be considered 
here.) 
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In a parallel development, the full perturbation ex- 
pansion for K(t) was derived from periodic-orbit theory 
B E El The connection between K(t), which is 
a quantum-mechanical object, and classical periodic or- 
bits follows from Gutzwillcr's trace formula, which ex- 
presses K (t) as a double sum over periodic orbits 7 and 

i nana, 

Kit) = [^A^e^-^/H^t- ^V^.( 6 ) 

Here ^4 7 , S* 7 , and T 7 are the stability amplitude, clas- 
sical action, and period of 7, respectively. While the 
leading contribution to K{t) comes from the diagonal 
terms 7 = 7' (up to time-reversal, if fi = 1) 01 > on? ~ 
diagonal contributions are responsible for 5K(t). Sieber 
and Richter and Heusler ei al. 0, 0, EH succeeded in 
classifying the relevant off-diagonal orbit pairs and calcu- 
lated their contribution to 5K(t) in the limit te/^h - > 0. 

Below, we show that periodic-orbit theory can also 
be used to calculate the Ehrenfest-time dependence of 
5K(t). Interestingly, while we confirm the field-theoretic 
result for the 0{t 2 ) term in 8K\{t) |5l21|. our result for 
the leading 0(t 3 ) contribution to SKzft) differs from that 
of Ref . yj 

V 2 

SK 2 (t) = (7) 

x [G(t - 2r E ) - 2Q(t - 3r E ) + Q{t - 4r E )] . 

In particular, we find that the minimum duration of off- 
diagonal pairs of orbits that contribute to SK2 is 2te, not 
3te. We also find that there are no TE-dependent cor- 
rections for t > 4te, in contrast to Ref. 0, where the te- 
dependent corrections to 8K2 persist up to t ~ tfj- (How- 
ever, our leading-order perturbative calculation does not 
answer the question whether such TE-insensitivity at long 
times persist to the higher order terms in t/tn.) 

Instead of calculating K(t) directly, it is more conve- 
nient to calculate the Laplace transform 

K(a) = /^A 7 A;,e l ^- s v)A-«(^+^)/2\. (8) 

\7>7' / 

The leading diagonal contribution to K can be calculated 
using the sum rule of Hannay and Ozorio de Almeida [2^1 , 

so that 

K(a) = -^+5K(a). (10) 

The inverse Laplace transform of the first term in Eq. 
(|10|l reproduces the leading term in Eq. J2J) above. 
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FIG. 1: Schematic drawing of a pair of orbits (shown solid 
and dashed) contributing to the leading interference correc- 
tion to the spectral form factor in the presence of time-reversal 
symmetry (a) and in the absence of time-reversal symmetry 
(b,c,d). The true orbits are piecewise straight, with specular 
reflection off the cavity's boundaries. The small-angle self- 
encounters are shown thick. Panels (a) and (b) also show the 
definitions of various durations used in the text. 



The leading 0(t 2 ) quantum correction for K(a) ex- 
ists in the presence of time-reversal symmetry only. The 
relevant pairs of periodic orbits 7 and 7' are shown in 
Fig. The existence of such pairs was pointed out by 
Sieber and Richter 18]; an equivalent configuration of 
classical trajectories appears in the field-theoretic formu- 
lation |3, Il2l |21| and in the diagrammatic calculation of 
the form factor for disordered cavities Q. The periodic 
orbit 7 in Fig^has a small- angle self-intersection. There 
are two loops of duration T\ and T2 through which 7 re- 
turns to the self-intersection. The trajectory 7' is equal 
to 7 in one of these loops, whereas 7' is the time-reversed 
of 7 in the other loop. Following Ref. IE 0, 13 we per- 
form the sum over such periodic orbits with the help of 
a Poincare surface of section taken at an arbitrary point 
during the self-intersection. The Poincare surface of sec- 
tion is parameterized using stable and unstable phase 
space coordinates s and u, normalized such that dsdu is 
the cross-sectional area element. Denoting the coordi- 
nate differences between the two points where 7 pierces 
the Poincare surface of section by s and u, the action dif- 
ference S 1 — S 1 < = su [23l |3 . The duration t enc of the 
self-encounter is defined as the time during which the two 
stretches of 7 are within a phase space distance c, where c 
is a classical scale below which the classical dynamics can 
be linearized. The periodic-orbit sum is then expressed 
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in terms of an integral over s, u, T\ and T 2 

SK^a) = I dT l( lT 2 f dsdu^±^ 2 + 2tcnc)2 

J J-c ( 27r? 



x cos( S u/h)e~ a ^+™^ 

d 2 1 [ c , , cos(su/?i)e- 2Qt — , , 
dsdu — ro ^2, 7 ' ( n ) 



where 



W = (1/A)ln(c 2 /M)> 



(12) 



A being the Lyapunov exponent of the classical dynamics 
in the cavity. The factor i e nc in the denominator cancels 
an unwanted contribution from the freedom to choose the 
Poincare surface of section anywhere inside the encounter 
region. Whereas Refs. |y, 0, UJ3 calculate the remaining 
integral in the limit ar E — > 0, we need to consider the 
effect of a finite Ehrenfest time. Such integrals have been 
considered in the context of quantum transport and we 
will be able to obtain the remaining integral in Eq. 
as well as all other necessary integrals from the literature. 
The integral needed here has been calculated in Ref. |25| 
and gives 
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A ln 7f 
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This result is in agreement with the te dependence of SK\ 
calculated by Tian and Larkin Its inverse Laplace 
transform is Eq. (0} above. 

We now consider the Ehrenfest-time dependent correc- 
tion 5K 2 (t) to the spectral form factor in the absence of 
time-reversal symmetry. There are three classes of pe- 
riodic orbits that give a contribution to SK^if) to order 
t 3 . These are shown in Fig. [TJa-d. Now, the classical 
orbit 7 has two small-angle self encounters. Its partner 
orbit 7' follows 7 between the encounters, but connects 
the ends of the encounters in a different way. Figure 
^> shows two separate encounters. Figure ^ shows a 
"three-encounter" , which arises when the two small-angle 
encounters of Fig. QJd are merged along one of the con- 
necting stretches (while keeping the duration of the other 
stretches finite). The orbit 7 of Fig.^J; goes through three 
loops between returns to the encounter region. Finally, 
Fig.^i shows a periodic orbit for which the self-encounter 
fully extends along one of these loops |2(j. Note that 
the encounter region can not simultaneously extend along 
two or more of the loops, because then these loops, and, 
hence, the trajectories 7 and 7' would be equal. 

The configurations of Fig. [Jd and c were also consid- 
ered by Heusler et al. 0, 0, Ujj ■ These two contributions 
cancel in the limit oite — > 0, so that one finds 6K2 = 
in that limit. The contribution of the trajectories shown 



Fig. vanishes in the limit arp, — > 0, which is why it 
was not considered in Refs.l6ll7l.ll3. However, as we will 
show below, it is needed to be taken into account when 
calculating TE-dependent corrections. The field-theoretic 
calculation of Ref. also has two contributions to 5K 2 
only, one with two two-encounters and the other with a 
single three-encounter. 

The contribution SK 2 i, of the two separate encounters 
for the periodic-orbit pair in Fig.^ 3 factorizes. Taking a 
Poincare surface of section at each of the encounters and 
proceeding as in the calculation of 5K\, one finds 



5K 2 b(a) 
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where we included a combinatorial factor 1/4 to account 
for permutations of the 2x2 passages of 7 through the 
two encounters 0, El- This result is in agreement 
with that obtained by calculating the two two-encounter 
diagram in the field-theoretic approach [27j • 

For the calculation of the contribution of the trajectory 
of Fig.^; one needs only one Poincare surface of section, 
taken at a point where all three stretches of 7 are within 
a phase space distance c. Labeling the phase space co- 
ordinates of the three piercings of 7 through the surface 
of section as (si,Uj), i = 1,2,3, the action difference is 



S-y S-y 



(16) 



where s — si — S3, s' = Si — S5, u = u\ — 113, and 
u' = ui — 115. Integrating over the durations of the three 
stretches of 7 that connect the three-encounter to itself, 
we find 
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dsds dudu 



su + s u 



)/h], 



{2nh)H 2 i t' 



where t cnc is the duration of the encounter, 



t — — In 



,(17) 



A min(|s|, \s'\, \s — s'\) min(|u|, \u + u'\) 

and t' cnc is the time that all three trajectories involved in 
the encounter are within a phase space distance c, 



A max(|s|, |s'|, \s — s'\) max(|u|, \u'\, \u + u'\) 



The prefactor 1/3 in Eq. I|17fl accounts for permutations 
of the three passages of 7 through the three encounter. 
The integration domain in Eq. Ijl7|l is max(|u|, \u + 
u' |), max(|s|, |s'|, \s — s'\) < c. The exponential factor 
contains the total time 2t enc +t' enc that the orbit 7 spends 



4 



in the encounter region. Taking the remaining integral 
over s, s', u, and v! from Ref. l28l we find 



SK 2c (a) 
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d 2 1 



2wM 2 i da 2 a 2 



3e 



-3qte 
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-4aT B l 



.(18) 



Although this result agrees with the field-theoretic cal- 
culation in the limit cute — ► 0, 0| , the appearance of a 
term proportional to exp(— 4«te) -oscillation differs from 
the field-theoretic approach p|, in which only one term 
proportional to exp(— 3q;te) appears with prefactor —1 

Finally, we have to calculate the contribution from tra- 
jectories with a three-encounter where the encounter re- 
gion fully wraps around one of the loops. In order to 
make optimal use of the literature on the Ehrenfest-time 
dependence of quantum transport, we calculate this con- 
tribution in an indirect way: We again consider the case 
of two two-encounters, but now allow the encounters to 
approach each other and overlap along two pre-assigned 
stretches of 7. This situation shown in Fig. \T]p (again), 
and we allow the encounters to approach each other along 
the central loop in the figure. We take a Poincare surface 
of section at each encounter, and measure the durations 
of between the two surfaces of section along the central 
loop by T[ and T 2 . Along the remaining stretches (the 
outer loop in Fig. [TJj), we still require non-overlapping 
encounters in order to enforce 7^7'. Hence, we param- 
eterize their duration using times T3 and T4 measured be- 
tween the ends of the encounters. Finally, i en c,i and £ e nc,2 
denote time that the inner and outer loops are within a 
phase space distance c, and t s and t u are the durations 
of eventual stretches that the two segments of the outer 
loop are within a phase space distance c from each other 
but not from the inner loop |25|. (The times t s and t u 
are zero except in the case of overlapping encounters.) 
With this parameterization, the total duration of 7 is 

T 7 = T{ + T' 2 +T 3 +T 4 + t cncA + < cnCi2 + 2t s + 2t u . (19) 

We now consider the integral 



dT[dT 2 dT 3 dT 4 / dsidu\ds 2 du 2 



T^cos[(uiSi 



u 2 s 2 )/K\t 
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where Si and m are phase space coordinates at the 
two Poincare surfaces of section, i = 1,2. This in- 
tegral contains both the case that the two encounters 
are separate and the case that the two encounters over- 
lap. If the two encounters are separate, which requires 
T[ > A-Mnduisal/c 2 ) and T 2 > A" 1 ln(|u 2 si |/c 2 ), the 
two stretches in the outer loop are never close to each 
other, hence t s = t u = 0. One then recovers the expres- 
sion for 5K 2 f )1 multiplied by four because of the combi- 
natorial factor 1/4 which is present in Eq. i|15fl but not 
in Eq. Q2U[I. If the two encounters overlap at one end 
but not at the other end, one recovers the scenario for 



SK 2c , multiplied by three because of the combinatorial 
factor 1/3 which is present in Eq. I|17|) but not in Eq. 
(I20JI . (Note that in this case the times t s and t u need not 
be zero.) If the encounters overlap at two ends, they span 
the reference loop. This scenario is neither contained in 
5K 2 b nor in 5K 2c . Since there is no combinatorial factor 
in this case, this is precisely the contribution 5K 2( i cor- 
responding to the trajectories of the type shown in Fig. 
^i. Hence 



SK 2d {a) =1- UK 2b (a) - 36K 2c (a). 



(21) 



From Sec. IV of Ref. |2fj, where an integral similar to / 
was calculated, we find 



Combining everything, we arrive at 



(22) 



5K 2 (a) 



5K 2b {a) + 6K 2c (a) + SK 2d {a) 
3 d 2 e - 2ar * 



2itU\ da 2 



[1 



(23) 



The inverse Laplace transform of this result is Eq. (JJJ 
above. 

The main difference between our result and that previ- 
ously obtained by the field-theoretic approach is that, in 
contrast to Ref. |5|, in Eq. (JJJ univers al q uantum correc- 
tions already appear after a time 2te 29] . This shortest- 
duration contribution to SK 2 stems from periodic orbit 
pairs which (for a part of their duration) wind around 
another, shorter, periodic orbit. Such orbits explained 
[25| the numerically observed TE-independence of con- 
ductance fluctuations in chaotic cavities H3, E3| An- 
other difference between our result and that previously 
obtained by the field-theoretic approach is the appear- 
ance of a discontinuity at 4te for the contribution SK 2c 
of a single three encounter. A feature at 4te, which is 
absent in the field-theoretic calculation, is essential for 
the validity of the 'effective random matrix theory' ^{| 
of the Ehrenfest-time dependence of the spectral gap in 
a chaotic cavity coupled to a superconductor psl ]. 

The calculation of the TE-dependence of the pertur- 
bative correction 8K 2 is, to the best of our knowledge, 
the first example in which the periodic orbit and the 
field-theoretical approaches lead to different results in the 
universal range. It is also the first example in which the 
two approaches are compared for a quantity that involves 
more than a single two-encounter. We believe this dis- 
crepancy can be traced to a difference in the way both 
approaches handle the encounter regions, in particular to 
differences implicitly or explicitly imposed by the short- 
time regularization procedures. 

In the periodic-orbit approach, a generic phase-space 
density function is described in reference to a classical 
trajectory. It expands and contracts along the unsta- 
ble and stable directions in phase space, but with a 
uniform density function at each point during the en- 
counter. The total duration of the encounter region is 
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~ te, beyond which the density functions are projected 
onto the low-lying sectors of the Perron- Frobenius modes 
and evolve irreversibly , despite the deterministic na- 
ture of the classical dynamics. On the other hand, in the 
field-theoretic approach of Ref. |5j, each 'orbit' represents 
the evolution from a singular density field concentrated 
at a point in phase space. Under the action of ellip- 
tical regulators, such singular fields first become locally 
smooth, spreading over the scale determined by the regu- 
lator strength. This scale is adjusted to be ~ Ti in order to 
mimic the evolution of quantum wave packets with mini- 
mal variance |H Ull Il4| | . At the time te/2 away from the 
origin (both forward and backward in time), the (now) 
locally smooth fields are projected onto the low-lying sec- 
tors of the Perron- Frobenius modes (2*3 . The two ap- 
proaches give the same result for a single two-encounter: 
in both cases the encounter region is a generic classical 
trajectory of duration rg. However, whether encounters 
are built from a reference trajectory or a reference 'point' 
matters for encounters of more than two trajectories. In 
the periodic-orbit approach, a reference orbit can be a 
periodic orbit with an arbitrarily short period. These 
reference orbits then give the universal contribution to 
SK at the time 2te that was calculated above. The ref- 
erence orbit can also be a trajectory of length < 2te, 



which is involved in two almost non-overlapping two en- 
counters. This case gives the feature at time 4te for the 
three-encounter contribution to the form factor. In the 
field-theoretic approach of Ref. 0, the evolution of the 
phase space density field is always altered irreversibly 
a time te/2 away from the reference point, thus, in the 
case of a three-encounter, leaving no room for a quantum 
correction at time 2te, nor for a feature at 4te- 

Summarizing, we used periodic orbit theory to calcu- 
late the leading Ehrenfest time dependent corrections 
to the spectral form factor K(t). We find that, in the 
absence of time-reversal symmetry, the corrections to 
the leading RMT prediction do not strictly vanish for 
2te < t < 4te if the Ehrenfest time te is finite. The 
fact that finite- te corrections give a nonzero correction 
to K(t) is in agreement with a previous field-theoretic 
calculation although our detailed expression for K(t) 
is not. We attribute the difference to the different short- 
time regularization used in Ref. 0- 
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